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Abstract — In a multicell multiuser MIMO downlink employing 
random beamforming as the transmission scheme, the heteroge- 
neous large scale channel effects of intercell and intracell interfer- 
ence complicate analysis of distributed scheduling based systems. 
In this paper, we extend the analysis in 1 1 1 and 1 2 1 to study the 
aforementioned challenging scenario. The cumulative distribution 
function (CDF)-based scheduling policy utilized in 1 1 1 and 1 2 1 is 
leveraged to maintain fairness among users and simultaneously 
obtain multiuser diversity gain. The closed form expression of 
the individual sum rate for each user is derived under the CDF- 
based scheduling policy. More importantly, with this distributed 
scheduling policy, we conduct asymptotic (in users) analysis to 
determine the limiting distribution of the signal-to-interference- 
plus-noise ratio, and establish the individual scaling laws for each 
user. 

I. Introduction 

With the emerging heterogeneous cellular structure [3J and 
the ever shrinking cell size, achieving high capacity with low 
design complexity in a multicell multiuser MIMO downlink 
has drawn considerable interest in recent years, e.g, see |S] 
and the references therein. Distributed scheduling policies are 
often favored due to operational scalability and affordable 
complexity incurred by the limited capacity of the backhaul. 
Under the employed distributed scheduling policy, analysis on 
the multicell network builds upon the extensive studies and 
insights drawn from the single cell network. 

For the single cell network without intercell interference, 
capacity boosting scheme relies on the independent varying 
channels across users, i.e., the well known multiuser diversity 
gain 111. To further harness this gain with multiple antennas, 
the notion of opportunistic beamforming is proposed J^, 
which is later extended to the notion of random beamforming 
Q to have same sum capacity growth as nonlinear precoding 
schemes IS), ||9| with reduced feedback requirement |10j. 
Multiuser diversity depends heavily on the scheduling policy, 
and it is important to guarantee scheduling fairness while 
achieving this gain in consideration of the heterogeneous 
large scale channel effects. This issue is tackled in [21 (the 
closed form sum rate in a homogeneous setup is derived 
in ifTTI '). by leveraging the cumulative distribution function 
(CDF)-based scheduling policy [12, 1 to satisfy the two desired 
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features: multiuser diversity and scheduling fairness. Accord- 
ing to the CDF-based scheduling policy |,12J . each user can 
be equivalently viewed as competing with other users with 
the same CDF, and thus each user's rate is independent of 
the statistics of other users. This interesting property enables 
a "micro" understanding of each user's rate performance 
compared to the conventional "macro" understanding of the 
sum rate performance. Due to this property, the notion of 
individual sum rate and individual scaling laws are proposed 
in |'2l to further understand both the exact and the asymptotic 
performance of random beamforming in a heterogeneous setup 
(and with different selective feedback schemes). 

In a multicell network, the heterogeneous channel effects 
come naturally from the different experienced intercell in- 
terference across users, even for a SISO setup. This issue 
is investigated in [11 for a generic single antenna based 
heterogeneous multicell OFDMA network, with both exact 
rate expression and asymptotic rate approximation derived. In 
ifTsl . the rate of convergence and the individual scaling laws 
are established for the SISO multicell setup. Both [JJ and 
ifTSl employ the CDF-based scheduling to maintain the two 
aforementioned scheduling features. In [14], the rate scaling 
for a power controlled network is examined with additional 
distance-based random variables. Since the time variations 
for the large scale and small scale channel effects are vastly 
different, both [11 and |13| concentrate on the randomness of 
the small scale channel effects. In |15||, a normalized form 
of transformation is applied in a multicell network to achieve 
fairness. The main difference of using CDF-based scheduling 
lfT2l is the inherent nonlinear functional transformation to 
strictly guarantee user fairness. 

In this work, we extend the analysis in ||T] and |l2| to a 
generic multicell multiuser MIMO setup. The random beam- 
forming is utilized as the multi-antenna transmission scheme 
to reduce feedback need (for part of the literature survey 
regarding random beamforming, please refer to [|2]). With 
spatial multiplexing in each cell, both intercell interference 
and intracell interference exist and users would experience 
heterogeneous interference. Under the CDF-based scheduling 
policy, we firstly derive the individual sum rate for each user 
from the exact analysis perspective. We further prove the type 
of convergence and the rate of convergence to the limiting 
distribution to establish the individual scaling laws. 



II. System Model 

Consider the downlink of a generic multicell multiuser 
MIMO network. We assume a narrowband model and the 
established analysis in this work can be extended to the wide- 
band model such as OFDMA using the techniques developed 
in f2\. Full spectrum reuse is assumed, and the process of cell 
association is assumed to be performed in advance. Without 
loss of generality, one base station Bq equipped with M 
antennas from the base station set B and its associated single 
antenna users ICq = {1, . . . ,k, . . . , Kq} with \ICo\ = Kq are 
considered. The random beamforming strategy at base station 
Bq employs M random orthonormal vectors ^J^-* G C^^^^ 
for m = 1,...,M, where the 0^°'' are drawn from an 
isotropic distribution independently every T (denoting the 
channel coherence interval for the block fading model) channel 
uses |7|. Denoting Sm\t) as the mth transmission symbol at 
time t, the transmitted vector s'^^\t) from base station _Bo at 
time t, is given as: 



scale channel gain G which may consist of path loss, antenna 
gain, and shadowing, varies usually on the order of seconds. 
Therefore, G is assumed to be known in advance by the 
system, through infrequent feedback; and the elements of h 
are modeled as complex Gaussian with zero mean and unit 
variance. Since for a given user k, the zj!^^^'s are identically 
distributed and correlated, the beam index m can be dropped 
in the expression of the CDF, which is derived in the following 
lemma. 

Lemma 1. The CDF of ^ can be expressed as 
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Proof: The main technique relies on the use of the 
moment-generating function (MGF) llT6l . Denote = 
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The received signal y^j^^ of user k (the time variable t is 



Etif^'EtU K'€'\'- Then the CDF of can be 
derived using the following procedure: 



dropped for notational convenience) is represented by 
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where the superscript indicates the base station with (0) being 
the cell of interest and (6 0) being the interfering cells. 
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denotes the number of interfering cells for user k, and v]^ where (a) follows from the fact that P Uf^ < x{ci°^ + 1] 

Hpnntps thf» aHHitivp whitp nnisp Histrihntprl with r'.\f(0 nr\ , . „ . ^ 



denotes the additive white noise distributed with CAf{0,af.). 
s^°) and s'''^ are the transmitted symbols by the serving 
cell and the interfering cell Bt, with E[|s(*'^|^] = po and 
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and h[^^ e Ci^*\ which are 



assumed to be independent across users, denote the small scale 
channel gains between the serving cell and user k, and between 
the interfering cell Bf, and user k. G^,"^ and g\^^ represent the 
large scale channel gains between the serving cell and user k, 
and between the interfering cell Bi, and user k respectively. 

Based on the aforementioned assumption, denoting as 
the SINR of user k for beam m, it can be expressed as: 



corresponds to the CDF of the exponential distribution at 
■^(Cfc*^^ + 1). From (|5]l, we note that the expression inside the 
integral corresponds to the MGF of denoted as ^',^(0) (t). 
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can be obtained below: 
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Now we examine the statistics of Z. 
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Note that the time 



variations for the large scale and small scale channel effects 
are vastly different ([T]. The variation of the small scale channel 
gain h occurs on the order of milliseconds; whereas the large 
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Combing Q and (|6]l yields F^m expressed in (|4|i. ■ 

The SINR will be fed baclH] and used for scheduling, which 
is pursued next. 

III. CDF-BASED Scheduling Policy 
AND Individual Sum Rate 

After receiving the SINR^°^ from user k for beam to, 
the scheduler is ready to perform scheduling. In order to 
guarantee scheduling fairness and obtain multiuser diversity 
gain, we employ the CDF-based scheduling policy [12J as 

'Full feedback wherein each user feeds back the SINR for M beams is 
assumed. The established results in this paper can be extended to different 
fornis of selective feedback using the techniques developed in |2|. 
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the opportunistic scheduling poUc}0. According to this policy, 
the scheduler will utilize F (o), and performs the following 

k 

functional transformation: 
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The transformed random variable z'f']^ is uniformly dis- 
tributed ranging from to 1, and can be regarded as the 
virtual received SINR of user k for beam m. The transformed 
random variables z]^\_Js are i.i.d. across users for a given 
beam, which enables the maximization at the scheduler side 
to perform scheduling in a fair manner Denoting fc*^ as the 
random variable representing the selected user for beam m, 
then: 

fcl = arg max Z^.^l . (8) 
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After user k*^ is selected per (|8]l, the scheduler utilizes the 
corresponding zf} for rate matching of the selected user 

Let xln' be the SINR of the selected user for beam to, and 
now consider the sum rate for a given base station Bq defined 
as follows: 
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From the aforementioned analysis, the sum rate can be for- 
mulated as: 
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where (a) follows from the sufficient small probability that 
multiple beams are assigned to the same user;(b) follows 
from the change of variable x = F (o) (t) and the fairness 

I'm '™ 

property of the CDF-based scheduling policy. 

The CDF-based scheduling enables a "micro" level under- 
standing of each user's rate performance, from both exact and 
asymptotic perspective. The individual sum rate for user k is 
defined in [21 as the the individual user rate multiplied by the 
number of users in cell Bq, namely: 
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For detailed motivation of the CDF-based scheduling as well as the 
rationale behind the notion of individual sum rate and individual scaling laws, 
please refer to 11] and |2|. 



Now employing the steps described below, the closed form 
expression for the individual sum rate is derived. 

Step 1 (PDF Decomposition): This step, utilized similarly in 
lH] and ||2|, is the essential step to decompose d{F^{o) (x))^" 

k 

into the following amenable input for Step 2. 
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Step 2 (Partial Fraction Expansion ^TT^J): This step manip- 
ulates part of ( fT2t for further integration. 
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where the expressions for and ■0^''^ are given on top of 
this page. 

Combining the outcomes of the two aforementioned steps, 

we can derive the closed form expression for rI^^ in the 
following procedure: 
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where Ji(a,/3,7) = /o°° [i+x)(p+x)-< '^^^ l2(a,/3,7) = 
/o°° (^+3^)7 dx. The calculation for Ii (a, /3, 7) and X2 (a, /3, 7) 
has been discussed in |1|, and their closed form expressions 
can be found in [1. (42)] and [1, (43)]. 

Up to now, we have performed exact analysis and derived 
the closed form results for the individual sum rate for an 
arbitrary selected user in a given base station. The derived 
results extend the exact analysis in [[IJ (multicell SISO setup) 
and 111 (single cell random beamforming), and can serve as 
a theoretical reference for evaluating the system performance 
under the CDF-based scheduling policy. In the next section, we 
will perform asymptotic analysis to evaluate the rate scaling 
laws for Rf\ which helps in understanding the asymptotic 
behavior of an individual user. 



Due to the complicated form of F (o) in dU, we need to find 
upper and lower bound for Wk-Kg by constructive methods, 
namely, w^^fo ^ Wk-.Ko — "iw^ /f ^ , by deriving upper and lower 
bound for F„(o). Define b'^'" = arg min pl^\ and b'^^^ = 

k O'^b'^Jk 

arg max p'l^K Then, one suitable upper bound and lower 
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in the following lemma. 
Lemma 3. 
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IV. Individual Scaling Laws 

This section is devoted to the asymptotic analysis. Section 
IIV-AI shows the type of convergence of In Section lEB 
the convergence rate to the limiting distribution is studied to 
establish the individual rate scaUng laws. 

A. Type of Convergence 

Firstly, we need to examine the tail behavior of the statistics 

of which has the form presented in (|4]l. Tools from 

extreme value theory [TSl are to be utilized. The following 

lemma describes the tail behavior of i^„(o) . 

^k 

Lemma 2. F„(o) belongs to the domain of attraction of the 
Gumbel distribution, i.e., F (q) e P^Gs). 

^k 

Proof: (Sketch) In order to prove that F (o) G ^'(Gs), 
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equivalent condition is: lim / w^*" = 1. Since 

similar methodologies in proving [1 Corollary 1] can be 
employed, we omit the detailed proof. ■ 

B. Rate of Convergence 

Knowing the type of convergence can lead to asymptotic 
approximation for the individual sum rate, which is investi- 
gated in JT]. Herein, dealing with higher order moments of 
the extreme order statistics and the rate of convergence is 
of interest. To establish the convergence rate to the limiting 
distribution for an individual user, the following definition of 
the so called growth function fTF| is needed: 
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One important step in proving the rate of convergence is to 
solve for a suitable coefficient sequence Wk-Ko by solving the 
following equation: 
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Proof: (Sketch) The upper and lower bound can be 
constructed by examining the large scale channel effects of 
intercell interference for user k. Herein, we only provide an 

intuitive explanation. F'-'fa) can be obtained by assuming that 

^k 

the intercell and intracell interference has the same large scale 
channel effects p^'''' \ The F^fo) can be found by using a 
similar line of argument. ■ 
Employing Lemma [3] the upper and lower bound for wu-Kq 
can be derived and are provided in the following corollary. 
Corollary L 

-wtX = pT logi^o - pf\{.h + l)M - 1) log (pi'™"' logi^o) 
+ O(logloglogifo). (19) 

^'-'ko = log - pt\{Jk + 1)M - 1) log (pI"™"' log Ao) 
+ 0(logIogIogA:o). (20) 

Proof: The w^'^Kq obtained via solving 1 — 

F^%-,{w^\^^) — -j^. Substituting the expression of F'-^o, and 

taking the log operator of both sides yields: 
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Since w)^\^^ ^ 00 as Kq — > cxd, dominates ( |2TI ) and 
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SO we have w^\^^ ^ p^f'' logiiTo- Then the sequence of w^\^ 
can be further written as w^.^ = pi°' log Kq + d^.^ , where 
4Koeo(logA'o). 

The expression for rf^.^o can be derived by substituting the 
formulation of w^.^^ into (|2T]) . which solves d^./j-^ as 
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Therefore, w^.^^ exhibits a scaling performance as ex 
pressed in ( fT9] l. The corresponding analysis for 
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be conducted following the same line of arguments. ■ 
Using the results from Corollary [T] and observing the fact 
that log pj^'' or log p^. *' are inconsequential when Ki^ goes 
large, we have the following expression for Wk-Ko'- 
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Once the expression of Wk:Ka is obtained, we can have the 
following inequality for the selected user's SINR for beam m 
(denoted by Xm^ in Section HIB by employing Q Corollary 
A.l], as follows: 
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Now we can state the following theorem for the rate scaling 
performance of i?^""*. 
Corollary 2. 
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< M¥ {x(? < Wk:Ko + loglogifo} 

X log (^l + Wk:K„ + Pk ^ loglogi^o 

+ MP > Wk:Ko + pT loglogi^o} log(l + pTk^) 

< Af log (l + Wk:Ko + Pk^ log log Xo) + 0(1). (26) 

A lower bound for i?^"' can be obtained similarly. Thus the 
individual scaling law for R^^'' exhibits a AI log log Kq growth 
in the large user regime. ■ 
Remark: Corollary |2] informs us that after we bring in 
the notion of individual scaling law corresponding to the 
individual sum rate for an arbitrary selected user, it exhibits 
a Af loglogA'o growth. This property is desirable from the 
perspective of opportunistic scheduling (e.g., greedy schedul- 
ing). In addition to this property, with the nonlinear functional 
transformation, the CDF-based scheduling policy can guaran- 
tee scheduling fairness in terms of long term user fairness (i.e., 
each user is equiprobable to be scheduled irrespective of their 
intercell and intracell interference). 

V. Conclusion 

The analytical impact of CDF-based scheduling policy in 
two special scenarios (multicell multiuser SISO and single 
cell multiuser MIMO) has been investigated in (T] and Q. 
This work extends and generalizes our previous works by 
addressing the rate performance in a generic multicell mul- 
tiuser MIMO downlink, with random beamforming as the 
signal transmission scheme. The most challenging part of 



this generaUzation lies in the existence of both intercell and 
intracell interference, as well as their accompanying heteroge- 
neous large scale channel effects. The CDF-based scheduling 
helps us to deal with this challenging scenario, and enables a 
"micro" understanding of the rate performance for any selected 
user The closed form individual sum rate is derived employing 
the MGF and the PDF decomposition. With the constructed 
bounding technique, we also establish the individual scaling 
laws to show that CDF-based scheduling exhibits the same 
scaling performance as opportunistic scheduling (but achieves 
scheduling fairness additionally). 
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